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Generalized Space and Frequency Index Modulation 

T. Datta, H. S. Eshwaraiah, and A. Chockalingam, 


Abstract —Unlike in conventional modulation where informa¬ 
tion bits are conveyed only through symbols from modula¬ 
tion alphabets defined in the complex plane (e.g., quadrature 
amplitude modulation (QAM), phase shift keying (PSK)), in 
index modulation (IM), additional information bits are conveyed 
through indices of certain transmit entities that get involved in 
the transmission. Transmit antennas in multi-antenna systems 
and subcarriers in multi-carrier systems are examples of such 
transmit entities that can be used to convey additional informa¬ 
tion bits through indexing. In this paper, we introduce generalized 
space and frequency index modulation, where the indices of active 
transmit antennas and subcarriers convey information bits. We 
first introduce index modulation in the spatial domain, referred 
to as generalized spatial index modulation (GSIM). For GSIM, 
where bits are indexed only in the spatial domain, we derive the 
expression for achievable rate as well as easy-to-compute upper 
and lower bounds on this rate. We show that the achievable 
rate in GSIM can be more than that in spatial multiplexing, 
and analytically establish the condition under which this can 
happen. It is noted that GSIM achieves this higher rate using 
fewer transmit radio frequency (RF) chains compared to spatial 
multiplexing. We also propose a Gibbs sampling based detection 
algorithm for GSIM and show that GSIM can achieve better 
bit error rate (BER) performance than spatial multiplexing. For 
generalized space-frequency index modulation (GSFIM), where 
bits are encoded through indexing in both active antennas as 
well as subcarriers, we derive the achievable rate expression. 
Numerical results show that GSFIM can achieve higher rates 
compared to conventional MIMO-OFDM. Also, BER results show 
the potential for GSEIM performing better than MIMO-OEDM. 

Index Terms —Multi-antenna systems, multi-carrier systems, 
spatial index modulation, space-frequency index modulation, 
achievable rate, transmit RE chains, detection. 


I. Introduction 

Multi-antenna wireless systems have become very popular 
due to their high spectral efficiencies and improved perfor¬ 
mance compared to single-antenna systems m- 0. Practical 
multi-antenna systems are faced with the problem of maintain¬ 
ing multiple radio frequency (RF) chains at the transmitter and 
receiver, and the associated RF hardware complexity, size, and 
cost ||4l. Spatial modulation, a transmission scheme which uses 
multiple transmit antennas but only one transmit RF chain, 
can alleviate the need for multiple transmit RF chains 0- 
0. In spatial modulation, at any given time, only one among 
the transmit antennas will be active and the other antennas 
remain silent. The index of the active transmit antenna will 
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also convey information bits, in addition to the information 
bits conveyed through the conventional modulation symbol 
(e.g., chosen from QAM/PSK alphabet) sent on the active 
antenna. An advantage of spatial modulation over conventional 
modulation is that, for a given spectral efficiency, conventional 
modulation requires a larger modulation alphabet size than 
spatial modulation, and this can lead to spatial modulation 
performing better than conventional modulation 0, ||9l. 

In this paper, we take the view that spatial modulation is an 
instance of the general idea of ‘index modulation’. Unlike in 
conventional modulation where information bits are conveyed 
only through symbols from modulation alphabets defined in 
the complex plane (e.g., QAM, PSK), in index modulation 
(IM), additional information bits are conveyed through indices 
of certain transmit entities that get involved in the transmis¬ 
sion. Transmit antennas in multi-antenna systems, subcarriers 
in multi-carrier systems, and precoders are examples of such 
transmit entities that can be used to convey information bits 
through indexing. Indexing in spatial domain (e.g., spatial 
modulation, and space shift keying which is a special case 
of spatial modulation) is a widely studied and reported index 
modulation technique; see Q and the references therein. Much 
fewer works have been reported in frequency and precoder 
index modulation techniques; e.g., subcarrier index modulation 
in ifTOl . ifm . ifT^ . lfT3l . and precoder index modulation in IITtII . 
The focus of this paper is twofold: i) generalization of the idea 
of spatial modulation, which we refer to as generalized spatial 
index modulation (GSIM), and ii) generalization of the idea of 
index modulation to both spatial domain (multiple-antennas) 
as well as frequency domain (subcarriers), which we refer to 
as generalized space-frequency index modulation (GSFIM). 

In spatial modulation, the choice of the transmit antenna 
to activate in a channel use is made based on a group of 
m bits, where the number of transmit antennas is rit = 2™. 
On the chosen antenna, a symbol from an M-ary modulation 
alphabet A (e.g., M-QAM) is sent. The remaining rit — 1 
antennas remain silent. Therefore, the achieved rate in spatial 
modulation, in bits per channel use (bpcu), is log 2 nt-|-log 2 M. 
The error performance of spatial modulation has been studied 
extensively, and it has been shown that spatial modulation can 
achieve performance gains compared to spatial multiplexing 
insi, Space shift keying is a special case of spatial mod¬ 
ulation El, where instead of sending an M-ary modulation 
symbol, a signal known to the receiver, say +1, is sent on the 
chosen antenna. So, the achieved rate in space shift keying 
is log 2 rit bpcu. In spatial modulation and space shift keying, 
the number of transmit RF chains is restricted one, and the 
number of transmit antennas is restricted to powers of two. 
The first contribution in this paper consists of generalization 
of spatial modulation which removes these restrictions nn- 
(H], an analysis of achievable rate, and proposal of a detection 
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algorithm. In generalized spatial index modulation (GSIM), 
the transmitter has nt transmit antenna elements and 
transmit RF chains, 1 < n^/ < rit, and out of n* antennas 
are activated at a time, thereby [log 2 (^‘ )J additional bits are 
conveyed through antenna indexing. Spatial modulation and 
spatial multiplexing turn out to be as special cases of GSIM for 
Urf = 1 and Urf = nt, respectively. We derive the expression 
for the achievable rate in GSIM and easy-to-compute upper 
and lower bounds on this rate. We show that the achievable 
rate in GSIM can be more than that in spatial multiplexing, 
and analytically establish the condition under which this can 
happen. It is noted that GSIM achieves this higher rate using 
fewer transmit RF chains compared to spatial multiplexing. 
We also propose a Gibbs sampling based detection algorithm 
for GSIM and show that GSIM can achieve better bit error 
rate (BER) performance than spatial multiplexing. 

In the second contribution in this paper, we introduce 
GSFIM which uses both spatial as well as frequency domain 
to encode bits through indexing. GSFIM can be viewed as a 
generalization of the GSIM scheme by exploiting indexing 
in the frequency domain as well. Index modulation that 
exploits the frequency domain alone - referred to as subcarrier 
index modulation (SIM) - has been studied in Qoi- m. 
These works have shown that OFDM with subcarrier index 
modulation (SIM-OFDM) achieves better performance than 
conventional OFDM, particularly at medium to high SNRs. 
These works have not exploited indexing in the spatial domain 
in MIMO systems. Our contribution addresses, for the first 
time, indexing both in space as well as frequency in MIMO 
systems. In particular, we (i) propose a signaling architecture 
for combined space and frequency indexing, (ii) study in detail 
its achieved rate in comparison with conventional MIMO- 
OFDM, and {in) show that better performance compared 
to that in conventional MIMO-OFDM can be achieved in 
the medium to high SNR regime. The proposed GSFIM 
system has N subcarriers, nt transmit antennas, and n^f 
transmit RF chains, 1 < Urf < nt- In the spatial domain, 
Urf out of nt transmit antennas are chosen for activation 
based on [log 2 (^* )J bits. In the frequency domain, in a 
space-frequency block of size n^/ x N, information bits are 
encoded in multiple sub-blocks where each sub-block is of size 
Urf xnf and ^ is the number of sub-blocks. We characterize 
the achievable rate in GSFIM as a function of the system 
parameters. We show that GSFIM can offer better rates and 
less transmit RF chains compared to those in conventional 
MIMO-OFDM. It is also shown that GSFIM can achieve better 
BER performance than MIMO OEDM. 

The rest of this paper is organized as follows. In Section HU 
we present the GSIM system model, and a detailed analysis 
of achievable rate and rate bounds in GSIM. We quantify rate 
gains and savings in transmit RE chains in GSIM compared 
to spatial multiplexing. The proposed detection algorithm 
for GSIM and its BER performance are also presented. In 
Section [Till we present the GSEIM system model, analysis of 
achievable rate in GSEIM, and BER performance of GSEIM. 
Conclusions and scope for future work are presented in Section 

El 



Fig. 1. GSIM transmitter. 

II. Generalized spatial index modulation 

In this section, we consider generalized spatial index mod¬ 
ulation (GSIM) which encodes bits through indexing in the 
spatial domain. In GSIM, the transmitter has nt transmit 
antennas and nrf transmit RE chains, 1 < n^/ < nt- In 
any given channel use, n^/ out of nt antennas are activated. 
Information bits are conveyed through both conventional mod¬ 
ulation symbols as well as the indices of the active antennas. 
Spatial multiplexing becomes a special case of GSIM with 
7irf = nt- We present an analysis of the achievable rates 
in GSIM, which shows that the maximum achievable rate in 
GSIM can be more than the rate in spatial multiplexing, and 
that too using fewer transmit RE chains. 

A. System model 

A GSIM transmitter is shown in Eig. [T] It has nt transmit 
antennas and n^/ transmit RE chains, 1 < n^-f < n*. An nrf x 
nt switch connects the RE chains to the transmit antennas. In a 
given channel use, n^/ out of nt transmit antennas are chosen 
and nrf M-ary modulation symbols are sent on these chosen 
antennas. The remaining nt — nrf antennas remain silent (i.e., 
they can be viewed as transmitting the value zero). Therefore, 
if A denotes the M-ary modulation alphabet used on the active 
antennas, the effective alphabet becomes Aq = A U 0. 

Define an antenna activation pattern to be a rij-length vector 
that indicates which antennas are active (denoted by a ‘1’ 
in the corresponding antenna index) and which antennas are 
silent (denoted by a ‘0’)- There are L = ) antenna 

activation patterns possible, and K = [log 2 bits are 

used to choose an activation pattern for a given channel use. 
Note that not all L activation patterns are needed, and any 2^ 
patterns out of them are adequate. Take any 2^ patterns out of 
L patterns and form a set called the ‘antenna activation pattern 
set’, §. Let us illustrate this using the following example. Let 
rij = 4 and nrf = 2. Then, L = (^ = Q, K = [log 2 6J = 2, 
and 2^ = 4. The six antenna activation patterns are given by 

{[ 1 , 1 , 0 , 0 ]^, [ 1 , 0 , 1 , 0 ]^, [ 0 , 1 , 0 , 1 ]^, [ 0 , 0 , 1 , 1 ]^, 

[ 0 , 1 , 1 , 0 ]^,[ 1 , 0 , 0 , 1 ]^}. 

Out of these six patterns, any 2^ = 4 patterns can be 
taken to form the set S. Accordingly, let us take the antenna 
activation pattern set as 

S = {[1,1,0,0]^, [1,0,1,0]^, [0,1,0,1]^, [0,0,1,1]^}. 

Table U shows the mapping of data bits to GSIM signals for 
nt = 4, nrf = 2 for the above activation pattern set. Suppose 
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Data bits 
K = 2 

Antenna activity 
pattern 

Antenna status 

Ant.l 

Ant.2 

Ant. 3 

Ant.4 

0 0 

[1, 1, 0, 0] '' 

e A 

G A 

OFF 

OFF 

0 1 

[1, 0, 1, 0]''' 

e A 

OFF 

G A 

OFF 

1 0 

[0, 1, 0, 1] '' 

OFF 

G A 

OFF 

G A 

1 1 

[0, 0, 1, 1] ''' 

OFF 

OFF 

G A 

G A 


TABLE I 

Data bits to GSIM signal mapping for nt = 4, n^f = 2. 

A: M- ARY MODULATION ALPHABET. 


4-QAM is used to send information on the active antennas. 
Let X G Aq* denote the n*-length transmit vector. Let 010011 
denote the information bit sequence. GSIM translates these bits 
to the transmit vector x as follows: i) the first two bits are used 
to choose the activity pattern, ii) the second two bits form a 
4-QAM symbol, and in) the third two bits form another 4- 
QAM symbol, so that, with Gray mapping, the transmit vector 
X becomes 

X = [1-fj, 0, -1-j, 0]"^, 
where j = a/— 1. 


B. Achievable rates in GSIM 


The transmit vector in a given channel use in GSIM is 
formed using i) antenna activation pattern selection bits, 
and ii) M-ary modulation bits. The number of activation 
pattern selection bits is [log 2 number of M-ary 

modulation bits is Urf log 2 M. Combining these two parts, 
the achievable rate in GSIM with rit transmit antennas, n^f 
transmit RF chains, and M-QAM is given by 



Antenna index bits 


4- rirf log 2 M 
^ ^ > 
modulation symbol bits 


bpcu. 


( 1 ) 


Let us examine the GSIM rate i?gsim in (Uli in some detail. In 
particular, let us examine how Rgsim varies as a function of 
its variables. Fig. |2] shows the variation of Rg^im as a function 
of Urf for different values of rit = 4, 8,12,16, 22, 32, and 4- 
QAM. The value of Urf in the x-axis is varied from from 0 
to rit- As mentioned before, n^/ = rit corresponds to spatial 
multiplexing. The i?gsim versus rirf plot for a given rit shows 
an interesting behavior, namely, for a given rit, there is an 
optimum n^/ that maximizes the achievable rate Rgsim- Let 
R™^ denote the maximum achievable rate, i.e.. 


max (2) 

® l<n,.f<nt 

In Fig. 13 it is interesting to see that R^^ does not necessarily 
occur at n^/ = rit, but at some rirf < rit- i?gsim can exceed 
the spatial multiplexing rate of rit log 2 M whenever the first 
term in O exceeds {rit — rirf) log 2 M. The following theorem 
formally establishes the condition under which the R^-^ will 
be more than the spatial multiplexing rate of nt log 2 M. 

Theorem 1: The maximum achievable rate in GSIM is 
strictly greater than the rate achieved in spatial multiplexing 
(i.e., > nt log 2 M) iff nt > 2M. 

Proof: Consider the two terms on the right-hand side (RHS) 
of the rate expression O- The first term (contribution due to 
antenna index bits) increases when Urf is increased from 0 
to and then decreases, i.e., it peaks at n^/ = L^J- The 



Fig. 2. Achievable rate in GSIM, Rgsim^ ^ function of n^f for different 
values of nt, and 4-QAM. 


second term (contribution due to modulation symbol bits), on 
the other hand, increases linearly with n^/- These two terms 
when added can cause a peak at some in the range LfJ< 
rirf < rit- Observe that, as we reduce nrf below nt, we gain 
rate from the first term but lose rate in the second term. The 
rate loss in the second term is log 2 M bpcu per RF chain 
reduced. Therefore, we can rewrite O as 


R„ 


nt log2 M + 



-{nt - nrf)\og 2 M. 


(3) 


Case 1: nt > 2M 

If nt > 2M, then [log 2 nt\ > log 2 M. By putting nrf = 
nt — 1 in (H, we get 


-^gsim — nt log2 M -f [log2 ntj - log2 M. (4) 

Therefore, in this case, the Rgsim in (Hli is more than nt log 2 M, 
i.e., GSIM with nrf = nt — 1 RF chains achieves more rate 
than spatial multiplexing. This implies > nt log 2 M, 

i.e., the maximum rate available in GSIM is more than the 
spatial multiplexing rate. Conversely, if nt < 2M, we show 
below that R^^ is not more than the spatial multiplexing 
rate. 

Case 2: nt < 2M 
If nt < 2M, 


logs rit < 1 + logs (5) 

From the properties of binomial coefficients, we have 

("*) = ( ) 

\nrfj \nt - rirfj 

^ nt{nt-l)---{nrf + l) ^ 

l.2.---{nt-nrf) 2 ”*-”’-/-!' ^ 

Hence, 


log 2 < [{nt - n-rf) log, nt -nt + Urt + IJ (7) 

\nrf J _ 

< [(nt - nrf){l + logs M) - nt -I- Urf + IJ (8) 

< {nt — Urf + 1) logs Af (9) 

< (nt — n,./) logs M. (10) 
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The inequality in (|7]) is obtained by taking logarithm in (|6l), 
and ® is obtained from 0 and (|5]l. Hence, using 0, we 
obtain i?gsim < n-t log 2 M, for 1 < < rit, and thus, for 

rit < 2M, < rit log 2 M. Combining the arguments in 

Cases 1 and 2, we get Theorem 1. □ 

From Fig. [1] the following interesting observations can be 
made: 

1) by choosing the optimum {nt,nrf) combination (i.e., 
using fewer RF chains than transmit antennas, Urf < 
rit), GSIM can achieve a higher rate than that of spatial 
multiplexing where n^/ = rit, and 

2) one can operate GSIM at the same rate as that of spatial 
multiplexing but with even fewer RF chains. 

For example, for rit = 32, the optimum that maximizes 
i?gsim is 24 and the corresponding maximum rate, 
is 71 bpcu. Compare this rate with 32 log 2 4 = 64 bpcu 
which is the rate achieved in spatial multiplexing. This is a 
11% gain in rate in GSIM compared to spatial multiplexing. 
Interestingly, this rate gain is achieved using lesser number 
of RF chains; 24 RF chains in GSIM versus 32 RF chains 
in spatial multiplexing. This is a 25% savings in transmit RF 
chains in GSIM compared to spatial multiplexing. Further, if 
GSIM were to achieve the spatial multiplexing rate of 64 bpcu 
in this case, then it can achieve it with even fewer RF chains, 
i.e., using just 18 RF chains which is a 43% savings in RF 
chains compared to spatial multiplexing. Table |II] gives the 
percentage gains in number of transmit RF chains at achieved 
rate R = R^^^ and R = rit log 2 M, and the percentage gains 
in rates achieved by GSIM compared to spatial multiplexing 
for nt = 16, 32 with BPSK, 4-QAM, 8-QAM, and 16-QAM. 


C. Bounds on achievable rates in GSIM 

We now proceed to obtain bounds on the achievable rate in 
GSIM. From O, we observe that 


i?gsim < l0g2 


nt\ 


Ur flint - nrf)l 


+ nrf\og2M, (11) 


and 


i?gsim > log2 


nt'. 


+ Urf l0g2M—l. (12) 


^Urflint - nrf)l, 

From the properties of the factorial operator ||22l, we have 




< n! < 


=^( 7 )’ 


Vne N. 


(13) 


Let us define the function /(nt, n^/, log 2 M) as 

/(nt, Urf, log 2 M) = nt log 2 nt - n^f log 2 nrf 

-{nt - nrf) \ 0 g 2 int - nrf) + n^/ log 2 M. (14) 
Substituting (00 in (fTTI) . using (O, and simplifying, we get 


6 Tlf 

Rg,im < log 2 + 0.51og2 - -r 

ZTT Urf [Tit — Urf) 

+fint,nrf,log2M). 

In a similar way, using (fT3] l in (fT2l) . we can write 


(15) 


R. 


> log: 




0.51og2 


nt 


gsim / ^>.-62 2 ^ n f(n,-n f) 

C- / Lrj f't ' ^rj ) 

^fint,nrf,\og2M) - 1. 


Let us rewrite ( fTSl l and ( fT6l l in the following way: 

Rgsim ^ fri^t 1 nrf) ~\~ f 2 int,nrf) c±, (12) 

and 

Rgsim ^ /l (^i 7 ^r/) “t“ f2int,nrf) C2, (13) 

where /i(nt,n,/) = 0.5log2 / 2 (nt,n,/) = 

/(nt, nrf, log 2 M), Cl = log 2 and C 2 = log 2 ^ - 1. For 
a fixed nt, the maximum value of fi{nt,nrf) in the range 
1 < nrf < nt — 1 is obtained at n^/ = 1 or n^/ = nt — 1, 
and the maximum value is 0.51og2 Hence, 

max{/i(nt,nr/)} = 0.51og2 . (19) 

Also, the term /i(nt,nr/) is minimized for nrf = 
the minimum value is 0.51og2 ^ = 1 — 0.51og2nt for even 
nt, and is 0.51og2 > 1 — 0.51og2nt for odd nt- 

Hence, 


min{/i(nt,nr/)} > 1 - 0.51og2nt. (20) 

Therefore, from (fTOl l. (|20] | and ( fTSl l, (fT4l) . we obtain the upper 
bound on i?gsim as 

Rgsim < /(nt,nr/,log 2 M)+ 0.51og2 

Tit -L 

+ log2^- (21) 

In a similar way, from (fTbl l and (fT4l) . we obtain the lower 
bound on i?gsim as 


-Rgsim > /(nt,nr/,log2M) - 0.51og2nt 

+ log2^-- (22) 

Since nt, nrf and M take finite positive integer values, and 
because of the floor operation in the first term on the RHS in 
0, we can rewrite the bounds in (l2Tli and (l22Ti as 


-Rgsim ^ 


/(nt,nr/,log 2 M) + 0.51og2 


nt 


nt - 1 


■log2 


27r 


(23) 


and 


-^gsim ^ 


/(nt, nrf, log 2 M) - 0.5 log 2 nt 


-log. 


2 e2 


(24) 


Note that the above bounds on Rgsim can be computed easily 
for any nt, nrf, without the need for the computation of 
factorials of large numbers in the actual rate expression in 
O. Further, noting that the optimum n^/ that maximizes 
/ 2 (nt,nr-/) is given by 


ntM 


Irf 


(16) 


M + 1 


(25) 
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7l/-ary 

alphabet 

Percentage saving in no. of Tx 
RF chains at ij = 

_ 

Percentage saving in no. of Tx RF 
chains at i? = nt log 2 M 

Percentage increase in rate at 

D _ jymax 

^g™ 

nt = 16 

nt = 32 

nt = 16 

nt = 32 

nt = 16 

nt = 32 

BPSK 

31.25 

40.63 

68.75 

71.88 

43.75 

46.88 

4-QAM 

18.75 

25 

37.5 

43.75 

9.385 

10.94 

8-QAM 

6.25 

12.5 

18.75 

21.88 

2.08 

3.13 

16-QAM 

6.25 

3.13 

6.25 

9.38 

0 

0.78 


TABLE II 

Percentage saving in transmit RF chains and percentage increase in rate in GSIM compared to spatial multiplexing eor n* = 16,32 

AND BPSK, 4-/8-/16-QAM. 


we obtain upper and lower bounds on by substituting 

n*j in (|25T i into (|2^ and d24li . respectively, as 


TDinax ^ 
-^gsim — 


rit log2(M + 1) + 0.5 log2 


n-t 


Tit - 1 


■log2 


27r 


(26) 


and 


TDrnax \ 
-^gsim — 


/ nt, 


M 


n-t 


, v^' 

' I 0 S 2 — 5 - 


M + 1 


Aog^M - 0.5log2 nt 


(27) 


These bounds on can be calculated for any given n-t and 

M directly, without exhaustive computation of the rate for all 
possible values of rirf. From (|2^ and dZTl i. we observe that 
as n-t —>■ 00, can be approximated by nt log 2 (M + 1). 

Note that a spatial multiplexing system which uses a zero- 
augmented alphabet Aq achieves the rate of ntlog 2 (M -|- 1), 
if all the symbols in Aq are equiprobable. 

In Fig. |3(a)| we plot the upper and lower bounds of R^sim 
computed using (l23t and d24li . respectively, along with exact 
7?gsim, for n-t = 16 and BPSK (M = 2). The number of 
RF chains, n^f, is varied from 1 to 15. It can be observed 
that the upper and lower bounds are tight (within 2 bpcu of 
the actual rate). In Fig. |3(b)| we plot the upper and lower 
bounds of R^-^ obtained from (l26l l and (l27l i. respectively, 
for different values of nt and M = 2,4- (i.e., BPSK, 4- 
QAM). The corresponding exact values are also plotted 

for comparison. It can be observed that the lower and upper 
bounds of R^^ are within 2 bpcu of the exact R'^-^- 


D. GSIM signal detection 

In this subsection, we consider detection of GSIM signals. 
Let H denote the x n* channel matrix, where is the num¬ 
ber of receive antennas. Assume rich scattering environment 
where the entries of H are modeled as circularly symmetric 
complex Gaussian with zero mean and unit variance. Let y 
denote the x 1-sized received vector, which is given by 

y = Hx + n, (28) 

where x is the n* x 1-sized transmit vector and n is the x 
1-sized additive white Gaussian noise vector at the receiver, 
whose ith element ~ tr^), Vz = 1, 2, • • • , n^. Let U 

denote the set of all possible transmit vectors, given by 

U = {x|x e Ao"'‘'^\||x||o = e §}, (29) 




(b) 

Fig. 3. (a) Bounds on i?asim BPSK for nt = 16 and varying (b) 

Bounds on BPSK and 4-QAM for varying nt. 

where ||x||o denotes the zero norm of vector x (i.e., number of 
non-zero entries in x), and denotes the antenna activation 
pattern vector corresponding to x, where = l,iff Xj ^ 
0,Vj = 1,2, ••• ,nt. Note that |U| = 2'^s*™. The activation 
pattern set S and the mapping between elements of § and 
antenna selection bits are known at both transmitter and 
receiver. Hence, from (|28]) and ( l29b . the ML decision rule 
for GSIM signal detection is given by 

X = argmin ||y — Hx|p. (30) 

For small values of nt and n-rf, the set U may be fully 
enumerated and ML detection as per (l30l) can be done. But 
for medium and large values of nt and n^-f, brute force 
computation of x in ( |30] | becomes computationally prohibitive. 
Here, we propose a low complexity algorithm for detection of 
GSIM signals. 
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The proposed approach is based on Gibbs sampling, where 
a Markov chain is formed with all possible transmitted vectors 
as states. As the total number of non-zero entries in the 
solution vector has to be equal to Urf, one can not sample 
each coordinate individually as is done in the case of Gibbs 
sampling based detection in conventional MIMO systems ll^ . 
To address this issue, we propose the following sampling 
approach; sample two coordinates at a time jointly, keeping 
other {rit — 2) coordinates fixed which contain (urf — 1) non¬ 
zero entries. 

1) Proposed modified Gibbs sampler: For any vector 
€ Aq* , ||o = Urf, where the t in the superscript 

of x^*^ refers to the iteration index in the algorithm. Let 
denote the locations of non-zero entries and 
jii J 2 , • • • ,jnt-nrf denote the locations of zero entries in x^*). 
We will sample and jointly, keeping other coordinates 
fixed, where I = 1,2, ■ ■ ■ , n^/ and fc = 1, 2, • • • , {rit — rirf). 
As any possible transmitted vector can have only Urf non¬ 
zero entries, the next possible state can only be 

any one of the following 2|A| candidate vectors denoted by 
{z“, tc = 1, 2, • • • , 2|A|}, which can be partitioned into two 
sets. In the first set corresponding to tu = 1, 2, • • • , |A|, we 
enlist the vectors which have the same activity pattern as x^*). 
Hence, = A“', 2 :“ = 0,Zq = x^q\q = 1,2,-■■ ,nt,q 
ihjk^w = 1,2, • • • , |A|. For w = |A| -f 1, |A| -f 2, • • • ,2|A|, 
we enlist the vectors whose activity pattern differs from that 
of x^*) in locations jk and ii. Hence, = 

0, Zg = x^q\q = 1,2, ••• ,nt,g ^ ityjk^w = 1,2, ••• , |A|. 
For w = |A| -f 1, |A| -f 2, • • • ,2|A|, we enlist the vectors 
whose activity pattern differs from that of x^*^ in locations 
jfc and ii. Hence, = 0,Zq = x^q\q = 

1,2, - •• ,nt,qf=ii,jkyw=\A\ + 1, |A| -f 2, • • • ,2|A|. 

To simplify the sampling process, we calculate the best 
vectors from the two sets corresponding to not swapping and 
swapping the zero and non-zero locations, and choose among 
these two vectors. Let denote the best vector from the 
first set corresponding to no swap. We set = x*^*) -f Ae^, 
and minimize ||y — Hx^'^jp over A. For this, we have 

||y-Hx^^f = ||y-H(xW+Ae,,)f 

= y^y - 25R (y^^x^*)) + x^^^^'Rx^*) 

-25R(Ay'"^^e,,) +25R(Ax(‘)""Re,,) +|ApR,,,,,, (31) 


where = y^H and R = H^H. Differentiating (ISTT i 

w.r.t A and equating it to zero, we get 


^opt 


Rt; A] 


H 


(32) 


where is the i^h column vector of R. We obtain x^^ = 
[x^*^ +Aoptei,]A, where [x]a denotes the element-wise quanti¬ 
zation of X to its nearest point in A. Similarly, we obtain x^, 
the best vector from the second set corresponding to swap. The 
next state is chosen between x'® and x^'® with proba¬ 

bility and respectively, where = (1 — q)p^ + |, 


„NS 


= 1 — p*, and 


r = 


1 + exp(- lly-HxS|F Jy-Hx^S||2 ^ ■ 


( 33 ) 


Here, q gives the probability of mixing between Gibbs sam¬ 
pling and sampling from uniform distribution. We use q = ^, 
because the simulation plots of BER as a function of q have 
shown that the best BER is achieved at around q = —. After 
sampling, the best vector obtained so far is updated. The above 
sampling process is repeated for all I and k. The algorithm 
is stopped after it meets the stopping criterion or reaches the 
maximum number of allowable iterations, and outputs the best 
vector in terms of ML cost obtained so far. 

2) Stopping and restart criterion: The following stopping 

criterion and restart criterion are employed in the algorithm. 
Let us denote the best vector so far as z. The stopping 
criterion works as follows; compute a metric 0s(z) = 

max (cmin, Cl exp(^(z))) , where (()(z) = 
is the normalized ML cost of z. If z has not changed for 
0s (z) iterations, then stop. This concludes one restart and z 
is declared as the output of this restart. Now, check whether 

belongs to S or not to check its validity. Several such runs, 
each starting from a different initial vector, are carried out till 
the best valid output obtained so far is reliable in terms of ML 
cost. Let us denote the best vector among restart outputs as s 
and the number of restarts that has given s as output as rg. 
We calculate another metric 0r(s) = [max (0, C 2 ^!'(s))J -|- 1 
and compare Vg with this. If Vg is equal to 0r-(s) or maximum 
number of restarts is reached, we terminate the algorithm. The 
listing of the proposed algorithm is given in Algorithm [T] 

3) Complexity: The complexity of the proposed Gibbs 
sampling based detector can be separated into three parts; 
i) computation of starting vectors, ii) computation of y^^ 
and R, and Hi) computations involved in the sampling and 
updating process. In our simulations, we use MMSE output as 
the starting vector for the first restart, and random starting vec¬ 
tors for the subsequent restarts. The MMSE output needs the 
computation of (H^H + cr'^lnt) H^y, whose complexity 
is 0{n^). Note that this operation includes the computations 
of y^^ and R. Eor the sampling and updating process, in 
each iteration, i.e., for each choice of I and k, the algorithm 
needs to compute x^*) and x^*) r^,., which requires 
0{nrf) computations. The rest of the computations are 0(1). 
The number of iterations before the algorithm terminates is 
found to be 0{nrf{nt — Urf)) by computer simulations. 
Thus, the total number of computations involved in Hi) is 
0{n'^f{nt — nrf)). Hence, the total complexity of the proposed 
algorithm for GSIM detection is 0{nf) + 0(n^j{nt — nrf)). 


E. BER performance results 

We now present the BER performance of GSIM. Eor 
systems with small rit, we present brute-force ML detection 
performance. Eor systems with large rit where brute-force ML 
detection is prohibitive, we present the performance using the 
proposed detection algorithm. We also compare the perfor¬ 
mance of GSIM with the performance of spatial multiplexing. 
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Algorithm 1 Proposed Gibbs sampling based algorithm for 
GSIM detection_ 

1 : input: y, H, nt,nr/; MAX-ITR: max. no. of iterations; MAX- 
RST; max. no. of restarts; 

2 : Compute = y^H and R = H^H; initialize r = 0, 

K = 10 “, q=—\ 

3: : ML cost fn; 0s(.) : stopping criterion fn; 0r(.) : restart 

criterion fn; 

4: while r < MAX-RST do 

5: : initial vector G = Urf, P = 

<;i(x*'°^); z = x^°^; t = 0 ; 

6: while t < MAX-ITR do 

7: for 1 = 1 to Urf do 

8 : for fc = 1 to nt — Urf do 

9: find ii and jk indices; 

10 : Compute Aopt from 1321 ; compute x^® = [x^*^ + 

-^opteijA; compute x^®; 

11 : Compute from ll3^; compute p® = (l — g)p® + |, 

pNS 

12: Choose between x® and x^® with probability 

p® & p^®; 

13: 7 = 0(x(*+i)); 

14: if (7 < P) then 

15: z = x*'‘'''^b /3 = 7 ; calculate 0s (z); 

16: end if 

17: f = f + 1 ; = P\ 

18: end for 

19: end for 

20: if 0s(z) < t then 

21: if == then 

22 : goto step 26 

23: end if 

24: end if 

25: end while 

26: r = r + 1 ; 

27: if G § then 

28: \i p < K then 

29: K = P\ Ts = 1; s = z; Compute 0r(s); 

30: end if 

31: if P == K then 

32: Cs = rs + 1; 

33: end if 

34: if Ts == 0r(s) then 

35: goto step 39 

36: end if 

37: end if 

38: end while 

39: output: s. s : output solution vector 


For notation purpose, a GSIM system with rit transmit anten¬ 
nas and Urf transmit RF chains is referred to as “{nt,nrf)- 
GSIM” system. Also, we use the term “(n*, np/)-SM” system 
to refer the spatial multiplexing system where nt = tirf. The 
following parameters are used in proposed detection algorithm: 
Cmin = 10np/(nt - Urf), Cl = 10np/(nt - rip/) log2 M, 
MAX-ITR = 8ntnrf{nt - nrf)s/M, MAX-RST= 20, C 2 = 
0.5(1 + log 2 M). Let denote the minimum number of 
RF chains in GSIM that achieves the same rate as in spatial 
multiplexing for a given rit and M. Let denote the number 
of RF chains that achieves for a given rit and M. 

In Fig. a we show the BER comparison between i) (4,2)- 
GSIM with 4-QAM, ii) (4,1)-GSIM with 16-QAM, and Hi) 
(2, 2)-SM with 8-QAM, using = 2. Note that in all the 
three systems, the modulation alphabets have been chosen such 



Fig. 4. BER compai'ison between (4, 2)-GSIM, (4,1)-GSIM, and (2, 2)-SM 
systems with 6 bpcu, Ur = 2, and brute-force ML detection. 



Fig. 5. BER comparison between MMSE detection, proposed detection, 
and brute-force ML detection in (4,3)-GSIM and (8,7)-GSIM systems with 
rir = nt, and 4-QAM. 


that the rate is the same 6 bpcu. Since the systems are small, 
brute-force ML detection is used. It can be seen that (4,2)- 
GSIM system performs better than (2,2)-SM system. That is, 
for the same rate of 6 bpcu and = 2, GSIM achieves 
better performance than spatial multiplexing by about 1 dB 
better performance at 0.01 uncoded BER. As we will see in 
Eigs.|6]and|7l this improvement increases to about 1.5 to 2 dB 
for 24 bpcu and 48 bpcu systems. It is noted that GSIM needs 
extra transmit antennas than spatial multiplexing to achieve 
this improvement. But the additional resources used in GSIM 
are not the transmit RE chains (which are expensive), but only 
the transmit antenna elements (which are not expensive). It can 
also be seen that even (4,1)-GSIM performs close to within 
0.5 dB of (2,2)-SM performance in medium to high SNRs. 
This shows that GSIM can save RE transmit chains without 
losing much performance compared to spatial multiplexing. 

Fig- |5] shows the BER performance of different detection 
schemes for GSIM. (4,3)-GSIM and (8,7)-GSIM with = nt 
and 4-QAM are considered. Note that the choice of Urf in both 
systems corresponds to Three detectors, namely, MMSE 
detector, proposed detector, and brute-force ML detector are 
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considered. It can be seen that MMSE detector yields very 
poor performance, but the proposed detector yields a perfor¬ 
mance which almost matches the ML detector performance. 
The proposed detector achieves this almost ML performance 
in just cubic complexity in nt, whereas ML detection has 
exponential complexity in n*. 

In Lig. |6] we compare the performance of three systems, 
each achieving 24 bpcu: i) (8,8)-SM with 8-QAM and ML 
detection using sphere decoder (SD), ii) (12,8)-GSIM with 
4-QAM and proposed detection, and in) (12,12)-SM system 
with 4-QAM using generalized sphere decoder (GSDl^ All 
the three systems use Ur = 8. Lig. |6] shows that the (12,8)- 
GSIM with proposed detection outperforms (8,8)-SM with SD 
employing same RL resources by about 2 dB in high SNR 
regime by using four extra transmit antennas. The performance 
of (12,8)-GSIM with proposed detection is very close to that of 
(12,12)-SM system with GSD which uses more RL resources 
to achieve the same rate. Also, the proposed detector has 
a much lower complexity than GSD which has exponential 
complexity in rit- 

Lig. |2] shows the BLR comparison between GSIM and SM 
using same RL resources for rirf = Ur = rirf to achieve 
48 bpcu. GSIM uses nt = 22 and 4-QAM, whereas (16,16)- 
SM scheme uses 8-QAM modulation alphabet to match the 
rate. Lor GSIM, the proposed detection is used. Lor SM, sphere 
decoding is used. It can be seen that, (22,16)-GSIM scheme 
outperforms (16,16)-SM scheme using same RL resources by 
about 2 dB in the medium to high SNR regime by using six 
extra transmit antennas. Also, the proposed detection has a 
much lower complexity than SD. 

In Ligs. |6] and |7] we also observe that at low SNRs 
the SM schemes have better BLR performance compared to 
the corresponding GSIM schemes. This can be explained as 
follows. Lirst, it can be observed that, to achieve the same 
rate, GSIM needs smaller-sized constellation compared to SM. 
Hence, GSIM will have a larger minimum distance among the 
constellation points than that in SM. Second, unlike in SM 
where there are no antenna index bits, the following two types 
of error events are observed in GSIM: i) the antenna activity 
pattern itself is decoded wrongly, and thus both the antenna 
index bits and modulation symbol bits are incorrectly decoded, 
and ii) the antenna activity pattern is decoded correctly, but 
the modulation symbol bits are wrongly decoded. At medium 
to high SNRs, the error event of the second type is more likely 
to occur and therefore this type of error events dominates the 
resulting performance. Coupled with this, a larger minimum 
distance among constellation points in GSIM than that in the 
corresponding SM makes GSIM to outperform SM in medium 
to high SNRs. But at low SNRs, the error event of the first type 
is more likely to occur and this error event type dominates the 
resulting performance. Since there are no antenna index bits in 
SM, error events of the first type do not occur in SM, leading 
to better performance for SM in the low SNR regime. 

'since Ur = 8, the (12,12)-SM system is an underdetermined system. 
Therefore, we have used the GSD in El which achieves ML detection in 
such underdetermined systems. GSD for spatial modulation has been reported 
in 1251 . 



Fig. 6. BER comparison among three systems achieving 24 bpcu: i) (8,8)- 
SM system with 8-QAM, ii) (12,8)-GSIM system with 4-QAM, and Hi) 
(12,12)-SM system with 4-QAM, Ur = 8. 



Fig. 7. BER comparison between GSIM and SM systems using same RF 
resources for ttr = rtr/ to achieve 48 bpcu. 

III. Generalized Space-Lrequency Index 
Modulation 

In this section, we propose a generalized space-frequency 
index modulation (GSLIM) scheme which encodes bits 
through indexing in both spatial as well as frequency domains. 
GSLIM can be viewed as a generalization of the GSIM scheme 
presented in the previous section by exploiting indexing in the 
frequency domain as well. In the proposed GSLIM scheme, 
information bits are mapped through antenna indexing in the 
spatial domain, frequency indexing in the frequency domain, 
and M-ary modulation. After mapping, the signal is modulated 
using OLDM and is transmitted through the selected antennas. 
We obtain the rate equation for the proposed GSLIM system 
and study its achievable rate, rate variation as a function 
of the parameters involved, and the rate gain compared to 
conventional MIMO-OLDM. 

A. System model 

The proposed GSLIM system uses nt transmit antennas, 
nrf transmit RL chains, 1 < n^/ < nt, N subcarriers. 
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Fig. 9. Frequency indexing in GSFIM. 
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Fig. 8. Block diagram of GSFIM transmitter and receiver. 
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and rir receive antennas. The channel between each transmit 
and receive antenna pair is assumed to be frequency-selective 
fading with L multipaths. The block diagrams of the GSFIM 
transmitter and receiver are shown in Fig. [ 8 ] At any given 
time, only n^/ transmit antennas are active and the remaining 
nt — Urf antennas remain silent. The GSFIM encoder takes 
[log 2 bits and maps to Urf out of n* transmit antennas 

(antenna index bits). It also takes additional bits to index sub¬ 
carriers (frequency index bits) and bits for M-ary modulation 
symbols on subcarriers. The frequency and antenna indexing 
mechanisms are detailed below. 

1) Frequency indexing: Consider a matrix B of size Urf x 
N whose entries belong to Aq, where Aq = A U 0 with 
A denoting an M-wry modulation alphabet. The frequency 
index bits and M-ary modulation bits are embedded in B 
as follows. The matrix B is divided into nt sub-matrices 
Bi, B 2 , • • • B„^, each of size Urf ^ nf, where n/ = ^ is 
the number subcarriers per sub-matrix (see Fig. |9]l. Let k, 
1 < k < UrfUf denote the number of non-zero elements in 
each sub-matrix, where each of the non-zero elements belong 
to A. This fc is a design parameter. Then, for each sub-matrix, 
there are If = possible ‘frequency activation patterns’. 

A frequency activation pattern for a given sub-matrix refers to 
a possible combination of zero and non-zero entries in that 
sub-matrix. Note that not all I / activation patterns are needed 
for frequency indexing. Any patterns out of them, where 
h = [log 2 (" 1 "')] , are adequate. Take any patterns out 
of If patterns and form a set called the ‘frequency activation 
pattern set’, denoted by S/. The frequency activation pattern 
for a given sub-matrix is then formed by choosing one among 
the patterns in the set S/ using kf bits. These kj bits are the 
frequency index bits for that sub-matrix. So, there are a total of 
Ubkf frequency index bits in the entire matrix B. In addition 
to these frequency index bits, krib log 2 M bits are carried as 
M-ary modulation bits in the non-zero entries of B. 

Example: Let us illustrate this using the following example. 
Let Urf =2, N = 16, Ub = 4, and k = 7. Then, n/ = ^ = 4, 
= (®) = 8, fc/ = [log 2 8J = 3, and 2^f = 8. In this 
example, If = 2^f = 8, i.e., all the 8 possible patterns are in 
the frequency activation pattern set, given by 


Suppose A is 4-QAM. Let [00101001111000110] denote the 
information bit sequence for sub-matrix Bi. The GSFIM 
encoder translates these bits to the sub-matrix Bi as follows: 
the first 3 bits are used to choose the frequency activity pattern 
(i.e., 001 chooses the activation pattern [ 1111 ] 
set E>f above), and the next 14 bits are mapped to seven 4- 
QAM symbols so that one 4-QAM symbol gets mapped to 
one active subcarrier. The sub-matrix Bi then becomes 


-1-j 0 -1+j 1-j 

1-j -1+j -1-j 1+j 


where j = %/M. Likewise, the sub-matrices B^, i = 2, 3,4 are 
formed. The full matrix B of size Urf x N is then formed as 


B = [Bi B 2 B 3 B 4 ]. 


Each row of the matrix B is of dimension 1 x N. There are 
Urf rows. Each V-length row vector in B is fed to the lEET 
block in the OEDM modulator to generate an 7V-length OEDM 
symbol. A total of n^/ such OEDM symbols, one for each 
row in B, are generated. These Urf OEDM symbols are then 
transmitted through Urf active transmit antennas in parallel. 
The choice of these Urf active transmit antennas among the 
rit available antennas is made through antenna indexing as 
described below. 

2) Antenna indexing: The selection of Urf out of nt 
antennas for transmission is made based on antenna index bits. 
The antenna index bits choose an ‘antenna activation pattern’, 
which tells which Urf antennas out of n* antennas are used for 
transmission. There are la = {n*f) activation patterns 

possible, and fca = [ log 2 (^‘) J bits are used to choose one 
among them. These ka bits are the antenna index bits. Note 
that not all la activation patterns are needed, and any 2 *^“ 
patterns out of them are adequate. Take any 2^“ patterns out 
of la patterns and form a set called the ‘antenna activation 
pattern set’, denoted by Sq. 

Example: Let us illustrate this using the following example. 
Let nt = 3, = 2. Then, la = ( 2 ) =8, ka= [log 2 = 

[log 2 3J = 1, and 2 ^“ = 2 . The possible antenna activation 
patterns are given by {[1,1, 0]^, [1,0,1]^, [0,1,1]^}. The set 
§a is formed by selecting any two patterns out of the above 
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three patterns. For example, can be 

§, = {[ 1 , 1 , 0 ]^,[ 1 , 0 , 1 ]^}. 

An rirf X rit switch connects the transmit RP chains to 
the transmit antennas. The chosen Urf out of n* transmit 
antennas transmit the MIMO-OFDM symbol constructed using 
the frequency index bits and M-ary modulation bits. The 
active transmit antennas can change from one MIMO-OFDM 
symbol to the other. 


B. Achievable rate, rate variation, and rate gain 

In GSFIM, the information bits are encoded using i) fre¬ 
quency indexing over each sub-matrix B^, i = - ,nb, 

ii) M-ary modulation symbols in each sub-matrix, and in) 
antenna indexing. The number of frequency indexing bits per 
sub-matrix is [log 2 The number of M-ary modu¬ 

lation bits in each sub-matrix is k log 2 M. The number of 


antenna indexing bits is 


_l 0 g 2 iZ e 

parts, the achievable rate in GSFIM with rit transmit antennas, 
Urf transmit RF chains, N subcarriers, rib sub-matrices, and 
M-ary modulation is given by 


R. 


sfim — 



+ 


Ra 

f krib log 2 M 
\N+L-1 


bpcu. 


(34) 


Rq 


Note that in a conventional MIMO-OFDM system, there is 
no contribution to the rate by antenna or frequency indexing, 
and the achieved rate is only through M-ary modulation 
symbols. Also, in MIMO-OFDM, M-ary modulation symbols 
are mounted on all N subcarriers on each of the rirf active 
transmit antennas. Therefore, the achieved rate in MIMO- 
OFDM (with no antenna and frequency indexing) for the same 
parameters as in GSFIM is given by 

f^mimo-ofdm — izF+T — l^'^rfNhg^M bpcu. (35) 

From 01 and OS), we can make the following observa¬ 
tions: 

• conventional MIMO-OFDM becomes a special case of 
GSFIM for rirf = R-t, rif = N (i.e., rib = !)■ 

• GSIM presented in Section |II] becomes a special case of 
GSFIM for N = rif = rib = 1, k = rirf. 

• for rirf < nt, Ra > 0, which is the additional rate con¬ 
tributed by antenna indexing. In this case, in 01 

can be more or less compared to i?mimo-ofdm depending on 
the choice of parameters. For example, the parameter k 
can take values in the range 1 to rirfUf. An instance 
where Rgsnm is more than i?mimo-ofdm happens when 
k = rirf rif, in which case Rp = 0 and Rq = i?mimo-ofdm- 
Therefore, Ra is the excess rate (rate gain) in GSFIM 
compared to MIMO-OFDM. Likewise, an instance where 
-Rgsfim is less than i?mimo-ofdm happens when fc = 1 , in 



k 

(a) M = 2 



k 


(b) M = 4 


Fig. 10. Rate Ri = Rj: + Rq as a function of k for different values of 
Nf = rirf rif. 


which case becomes which is less 

than R ■ „ oivpn hv »f log 2 (M~) 

Ilian -rimimo-ofdm glVcH Dy 

• the sum of rates Rp and Rq in 01 as a function of 
k reaches its maximum for a value of k in the range 
^ rirfjif j rirfUf, and so does the total rate i?gsfim- The 
maximum Rgshm will be more than or equal to i?mimo-ofdm- 
We now illustrate the above observations through numerical 
results. Define i?i = Rp Rq and Nf = rirfUf. In Fig. 
[ini we plot i?i as a function of k, for different values of 
Nf = 8,16,32, L = 4, and M = 2,4:. We observe that Ri 
reaches its maximum value for k between and Nf. Also, 
the maximum i?i increases as Nf increases because the Rp 
term in 01 increases with Nf. 

In Fig. [TT] we plot the maximum i?gsfirn as a function of 
nt for Urf = 8, N = 32, L = 4, and n/ = 1,2,4,8,16,32. 
f?mimo-ofdm IS also plotted for Comparison. We observe that for 
a given nf, the maximum i?gsfim increases with n* because 
of the increase in antenna index bits carried. For a given n* 
and nrf, the maximum i?gsfim increases with increase in ny 
because of increase in Vy and the associated increase in Rp. 
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From this figure, we can see that GSFIM can achieve a rate 
gain of up to 65% for M = 2 and up to 19% for M = 4, 
compared to MIMO-OFDM. In Fig. [12] we have plotted the 
percentage rate gain in GSFIM compared to MIMO-OFDM 
(i.e., difference between maximum i?gsfim and i?mimo-ofdm in 
percentage), as a function of n^/ for rit = 32, N = 32, L = 4, 
and n/ = 2,4,8,16, 32. As can be observed in Fig.[T2| GSFIM 
can achieve rate gains up to 65% for M = 2 and 20% for 
M = 4, compared to MIMO-OFDM. 

In Fig. [Tsj we plot the maximum iigsHm as a function of 
rirf for a given nt = 32, N = 32, L = 4 and n/ = 1,32. 
We can observe that for a given n/, the rate increases with 
Urf because of the increase in i?F- For a given rirf, the the 
maximum Rg^fim increases with increase in n/. In Fig. [T4| 
we have plotted bar graphs showing the percentage savings 
in transmit RF chains in GSFIM compared to MIMO-OFDM 
Tit = N = 2)2, L = A, and rif = 1,4, 32. It can be observed 
that this savings is high for small-sized modulation alphabets 
- e.g., the savings is up to 42% for M = 2 and 20% for 
M = 4. 

In Fig. [T5| we plot the maximum i?gsfim as a function of 
Uf for nt = 32, n^/ = 8, L = A, and N = 32. We can 
observe that the maximum i?gsfim increases for up to certain 
Uf and thereafter it saturates. This is because the maximum 
Ri saturates to a value for large Nj. 

C. GSFIM signal detection and performance 

In this subsection, we consider GSFIM signal detection 
and performance. Let H„ denote x nt channel matrix on 
subcarrier n. Let H® denote the x n^/ channel matrix 
corresponding to the chosen antennas. The superscript a in 
refers to the antenna activation pattern that tells which 
antennas are chosen. Let us denote the Ur x 1-sized received 
vector on subcarrier n as y„, which can be written as 

y„ = H®z„-fw„, n = 1,2,--- ,iV, (36) 

where z„ is the n^/ x 1 -sized transmitted vector on subcarrier 
n, and w„ is the x 1-sized additive white Gaussian noise 
vector at the receiver, w„ ~ CA/'(0, Consider the 

system model in (l3^ for the jth sub-matrix, given by 

yi = Hfz;+Wi, I = ■ ■ ,inj, (37) 

where ij = {i — l)n/ A- j. Write (iTTl i as 

y* = Gfz*-bw*, 1,2, ••• ,n&, (38) 

where 


yii 


^ii 

yi2 


2^2 



L J 


r Hf, 0 1 

Hf. 

Gf = 

0 , . 

The ML metric for a given antenna activation pattern a and 
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(a) M = 2 


n^=8, N=32, L=4, M=4 

rf 




, GSFIM. n,=32 

GSFIM. n,=16 

GSFIM. n,=8 

GSFIM, n,=4 


^ GSFIM. nj=2 

GSFIM. n^-1 


_MIMO-OFDM, n =n^8 
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Fig. 11. Maximum Rgafim ^ function of nt, for n^f = 8 and different 
values of Uf. 

vectors z*,* = 1, • • • ,ni, representing the frequency activation 
pattern and M-ary modulation bits is 

rib 

d(a,z\z2,... ,z"^) = ^||y*-Gfz*f. (39) 

i=l 

Let U denote the set of all possible iVj^-length transmit vectors 
corresponding to a sub-matrix. Then, U is given by 

U = {x|xe Ao'^^^M|x||o = A:,C^gS/}, (40) 

where denotes the frequency activity pattern corresponding 
to X, where = l,iffxj f 0 ,Vj = 1,2,--- ,Nf. The 
antenna activation and frequency activation pattern sets (§o, 
§/), and the antenna and frequency index bit maps are known 
at both transmitter and receiver. Therefore, from ( |39] | and (l40l) . 
the ML decision rule for GSFIM signal detection is given by 

(a, z^, z^, • • • , z"'’) = argmin d(a, z^, z^, • • • , zf4).) 

By inverse mapping, the antenna index bits are recovered 
from a and the frequency index bits are recovered from 
Z^, Z^, • • • , z"*”. 
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(a) M = 2 



(b) M = 4 


Fig. 12. Percentage rate gain in GSFIM compared to MIMO-OFDM as a 
function of n^f and Uf. 


In Figs. [T6l a) and (b), we show the BER performance 
of GSFIM in comparison with MIMO-OFDM under MF 
detection. In Fig. [Tbl al. the GSFIM system has rit = 3, 
Urf = 2, N = 8, Uf = 4, = 2,4, 4-QAM, and 

the achieved rate is i?gsfim = 3.1818 bpcu. The MIMO- 
OFDM has nt = Urf = 2, N = 8, rir = 2, A, 4-QAM, 
and the achieved rate is i?mimo-ofdm = 2.9091 bpcu. In Fig. 
[Tbl b). the GSFIM system has nt = 3, Urf = 2, N = 16, 
Uf = A, Ur = 2,3, L = 4, 4-QAM, and the achieved 
rate is iigsfim = 3.6316 bpcu. The MIMO-OFDM system has 
nt = nrf = 2, N = 16, n^ = 2,3, L = 4, 4-QAM, and 
the achieved rate is i?mimo-ofdm = 3.3684 bpcu. It is seen 
that in Figs. [Tbl a) and (b), GSFIM has higher rates than 
MIMO-OFDM. In terms of error performance, while MIMO- 
OFDM performs better at low SNRs, GSFIM performs better 
at moderate to high SNRs. This performance cross-over can be 
explained in the same way as explained in the case of GSIM in 
the previous section (Sec. lII-El Figs. |6] and |7]i; i.e., at moderate 
to high SNRs, errors in index bits are less likely and this makes 
GSFIM perform better; at low SNRs, index bits and hence 




Fig. 13. Maximum -Rgstim ^ function of Urf, for nt = N = 32, and 
Uf = 1,4,32. 

the associated modulation bits are more likely to be in error 
making MIMO-OFDM to perform better. Similar performance 
cross-overs have been reported in the literature for single¬ 
antenna OFDM with/without subcarrier indexing (e.g., mx 
where it has been shown that OFDM with subcarrier indexing 
outperforms classical OFDM without subcarrier indexing at 
moderate to high SNRs, whereas classical OFDM outperforms 
OFDM with subcarrier indexing at low SNRs. The plots in 
Figs. [161 a) and (b) essentially capture a similar phenomenon 
when there are index bits both frequency as well as spatial 
domains. 

IV. Conclusions 

We introduced index modulation where information bits 
are encoded in the indices of the active antennas (spatial 
domain) and subcarriers (frequency domain), in addition to 
conveying information bits through conventional modulation 
symbols. For generalized spatial index modulation (GSIM), 
where bits are indexed only in the spatial domain, we derived 
the expression for achievable rate as well as easy-to-compute 
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(a) M = 2 



Fig. 14. Percentage savings in number of transmit RF chains in GSFIM 
compared to MIMO-OFDM, for nt = N = 32, rif = 1, A, 32. 



Fig. 15. Maximum -Rgsflni ^ function of for fixed nt, 




Fig. 16. BER performance of GSFIM and MIMO-OFDM under ML 
detection, (a) GSFIM with nt = 3, nj.f = 2, N = S, Uf = A, Ur = 2,4, 
L = A, 4-QAM, 3.1818 bpcu, and MIMO-OFDM with nt = nj.f = 2, 
N = S, nr = 2,A, L = A, 4-QAM, 2.9091 bpcu. (b) GSFIM with nt = 3, 
rirf = 2, N = 16, nf = A, nr = 2,3, L = A, 4-QAM, 3.6316 bpcu, and 
MIMO-OFDM with m = Urf = IQ, N = 16, nr = 2,3, L = 4, 4-QAM, 
3.3684 bpcu. 

upper and lower bounds on this rate. We showed that the 
achievable rate in GSIM can be more than that in spatial 
multiplexing, and analytically established the condition under 
which this can happen. We also proposed a Gibbs sampling 
based detection algorithm for GSIM and showed that GSIM 
can achieve better BER performance than spatial multiplexing. 
GSIM achieved this better performance using fewer transmit 
RF chains compared to spatial multiplexing. For generalized 
space-frequency index modulation (GSFIM), where bits are 
encoded in the indices of both active antennas as well as sub¬ 
carriers, we derived the achievable rate expression. Numerical 
results showed that GSFIM can achieve higher rates compared 
to conventional MIMO-OFDM. Also, BER results using ME 
detection showed the potential for GSFIM performing better 
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than MIMO-OFDM at moderate high SNRs. Low complexity 
detection methods for GSFIM can be taken up for future 
extension to this work. 
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